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Abstract — We describe a non-cooperative interference 
alignment (IA) technique which allows an opportunistic 
multiple input multiple output (MIMO) link (secondary) 
to harmlessly coexist with another MIMO link (primary) 
in the same frequency band. Assuming perfect channel 
knowledge at the primary receiver and transmitter, capac- 
ity is achieved by transmiting along the spatial directions 
(SD) associated with the singular values of its channel 
matrix using a water-filling power allocation (PA) scheme. 
Often, power limitations lead the primary transmitter to 
leave some of its SD unused. Here, it is shown that the 
opportunistic link can transmit its own data if it is possible 
to align the interference produced on the primary link with 
such unused SDs. We provide both a processing scheme 
to perform IA and a PA scheme which maximizes the 
transmission rate of the opportunistic link. The asymptotes 
of the achievable transmission rates of the opportunistic 
link are obtained in the regime of large numbers of 
antennas. Using this result, it is shown that depending on 
the signal-to-noise ratio and the number of transmit and 
receive antennas of the primary and opportunistic links, 
both systems can achieve transmission rates of the same 
order. 



I. Introduction 

The concept of cognitive radio is well-known by now. 
The main idea is to let a class of radio devices, called 
secondary systems, opportunistically access certain por- 
tions of spectrum left unused by other radio devices, 
called primary systems, at a given time or geographical 
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area |2|. These pieces of unused spectrum, known as 
white-spaces, appear mainly when either transmissions 
in the primary network are sporadic, i.e., there are 
periods over which no transmission takes place, or there 
is no network infrastructure for the primary system in 
a given area, for instance, when there is no primary 
network coverage in a certain region. In the case of dense 
networks, a white-space might be a rare and short-lasting 
event. As a matter of fact, the idea of cognitive radio 
as presented in O (i.e., spectrum pooling), depends on 
the existence of such white-spaces Q. In the absence 
of those spectrum holes, secondary systems are unable 
to transmit without producing additional interference on 
the primary systems. One solution to this situation has 
been provided recently under the name of interference 
alignment (IA) |4). Basically, IA refers to the construc- 
tion of signals such that the resulting interference signal 
lies in a subspace orthogonal to the one spanned by the 
signal of interest at each receiver Q. The IA concept 
was introduced separetely and almost simultaneously 
by several authors O, Q, 0, d. Recently, I A has 
become an important tool to study the interference 
channel, namely its degrees of freedom (H, (H, (9J. The 
feasibility and implementation issues of IA regarding 
mainly the required channel state information (CSI) has 
been also extensively studied 03, EH, D2L lPT3l . 
In this paper we study an IA scheme named opportunistic 
IA (OIA) HI. The idea behind OIA can be briefly 
described as follows. The primary link is modeled by 
a single-user MIMO channel since it must operate free 
of any additional interference produced by secondary 
systems. Then, assuming perfect CSI at both transmitter 
and receiver ends, capacity is achieved by implementing 
a water-filling power allocation (PA) scheme lfl4l over 
the spatial directions associated with the singular values 
of its channel transfer matrix. Interestingly, even if the 
primary transmitters maximize their transmission rates, 
power limitations generally lead them to leave some of 
their spatial directions (SD) unused. The unused SD 
can therefore be reused by another system operating 
in the same frequency band. Indeed, an opportunistic 
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transmitter can send its own data to its respective receiver 
by processing its signal in such a way that the inter- 
ference produced on the primary link impairs only the 
unused SDs. Hence, these spatial resources can be very 
useful for a secondary system when the available spectral 
resources are fully exploited over a certain period in a 
geographical area. The idea of OIA, as described above, 
was first introduced in 12 considering a very restrictive 
scenario, e.g., both primary and secondary devices have 
the same number of antennas and same power budget. In 
this paper, we consider a more general framework where 
devices have different number of antennas, different 
power budgets and no conditions are impossed over the 
channel transfer matrices (In 12, full rank condition was 
impossed over certain matrices). 

The rest of this paper is structured as follows. First, the 
system model, which consists of an interference channel 
with MIMO links, is introduced in Sec. [Ill Then, our aim 
in Sec. [HI] is twofold. First, an analysis of the feasibility 
of the OIA scheme is provided. For this purpose, the 
existence of transmit opportunities (SD left unused by 
the primary system) is studied. The average number 
of transmit opportunities is expressed as a function 
of the number of antennas at both the primary and 
secondary terminals. Second, the proposed interference 
alignment technique and power allocation (PA) policy at 
the secondary transmitter are described. In Sec. IIV-BI 
tools from random matrix theory for large systems are 
used to analyze the achievable transmission rate of 
the opportunistic transmitter when no optimization is 
performed over its input covariance matrix. We illustrate 
our theoretical results by simulations in Sec. [Vj Therein, 
it is shown that our approach allows the secondary link 
to achieve transmission rates of the same order as those 
of the primary link. Finally, in Sec. [VI] we state our 
conclusions and provide possible extensions of this work. 

II. System Model 

Notations. In the sequel, matrices and vectors are 
respectively denoted by boldface upper case symbols and 
boldface lower case symbols. An N x K matrix with 
ones on its main diagonal and zeros on its off-diagonal 
entries is denoted by InxK, while the identity matrix of 
size iV is simply denoted by Ijv- An N x K matrix with 
zeros in all its entries (null matrix) is denoted by OnxK- 
Matrices X T and are the transpose and Hermitian 
transpose of matrix X, respectively. The determinant of 
matrix X is denoted by |X|. The expectation operator 
is denoted by E[.]. The indicator function associated 
with a given set A is denoted by and defined 

by 1.4 (x) = 1 (resp. 0) if x € A (resp. x ^ A). 
The Heaviside step function and the Dirac delta function 



are respectively denoted by /i(-) and £(•). The symbols 
M, R, and C denote the sets of non-negative integers, 
real numbers, and complex numbers, respectively. The 
subsets [0, +oo[ and ] — oo,0] are denoted by R + and 
R~, respectively. The operator (x) + with x G II is 
equivalent to the operation max(0, x). Let A be an 
nxn square matrix with real eigenvalues \a,i, ■ ■ ■ , \a,u- 
We define the empirical eigenvalue distribution of A 
b y f a (') — ^Z)r=lM(<^ - ^A,i), and, when it exists, 
we denote (A) the associated eigenvalue probability 
density function, where Fa{-) and /a(0 are respectively 
the associated limiting eigenvalue distribution and prob- 
ability density function when n — > +oo. 
We consider two unidirectional links simultaneously op- 
erating in the same frequency band and producing mutual 
interference as shown in Fig. [2 The first transmitter- 
receiver pair (Txi,Rxi) is the primary link. The pair 
(Tx2,Rx2) is an opportunistic link subject to the strict 
constraint that the primary link must transmit at a rate 
equivalent to its single-user capacity. Denote by Ni 
and Mi, with i = 1 (resp. i = 2), the number of 
antennas at the primary (resp. secondary) receiver and 
transmitter, respectively. Each transmitter sends indepen- 
dent messages only to its respective receiver and no 
cooperation between them is allowed, i.e., there is no 
message exchange between transmitters. This scenario is 
known as the MIMO interference channel (IC) [Q3]], 
with private messages. A private message is a message 
from a given source to a given destination: only one 
destination node is able to decode it. Indeed, we do not 
consider the case of common messages which would be 
generated by a given source in order to be decoded by 
several destination nodes. 

In this paper, we assume the channel transfer matrices 
between different nodes to be fixed over the whole 
duration of the transmission. The channel transfer matrix 
from transmitter j G {1,2} to receiver i £ {1,2} is 
an Ni x Mj matrix denoted by Hy which corresponds 
to the realization of a random matrix with independent 
and identically distributed (i.i.d.) complex Gaussian cir- 
cularly symmetric entries with zero mean and variance 
■jjj:, which implies 

V(i, j) e {1,2} 2 , Trace (E [Hy H#]) = N { . (1) 

The Li symbols transmitter i is able to simultaneously 
transmit, denoted by s^i, . . . , s^z,., are represented by 
the vector s, = . . . , s^ij . We assume that 

Vi G {1,2} symbols s^i, . . . , are i.i.d. zero-mean 
circularly-symmetric complex Gaussian variables. In our 
model, transmitter i processes its symbols using a matrix 
Vj to construct its transmitted signal VjS,. Therefore, 
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the matrix Vj is called pre-processing matrix. Following 
a matrix notation, the primary and secondary received 
signals, represented by the JVj X 1 column-vectors y,- l , 
with i G {1,2}, can be written as 

/ n \ / Hn H 12 \ / Visi \,(ni\ 
\ r 2 J V H 21 H 22 J \ V 2 s 2 y + V n 2 ^ ' 

(2) 

where rij is an Ni -dimensional vector representing noise 
effects at receiver i G {1,2} with entries modeled by 
an additive white Gaussian noise (AWGN) process with 
zero mean and variance af, i.e.,Vi G {1, 2}, E [njnp] = 
ofl/v^ At transmitter i G {1,2}, the L, x Li power 
allocation matrix Pj is defined by the input covariance 
matrix Pj = E [sjS^] . Note that symbols s^i . . . , Si,z, ( , 
Vi G {1,2} are mutually independent and zero-mean, 
thus, the PA matrices can be written as diagonal matrices, 
i.e., Pj = diag (pi,i, . . . ,Pi,L t )- Choosing P.; therefore 
means selecting a given PA policy. The power constraints 
on the transmitted signals VjSj can be written as 

Vi G {1, 2} , Trace (V 4 P,Vf ) < M< Pi , max . (3) 

Note that assuming that the i.i.d. entries of matrices Hy, 
for all G {1,2} 2 , are Gaussian random variables 

with zero mean and variance j^-, together with the power 
constraints in ©, is equivalent to considering a system 
where the entries of matrices H^- for all G {1,2} 2 
are Gaussian random variables with zero mean and unit 
variance, and the transmitted signal VjSj are constrained 
by a finite transmit power Pi >ra ax- Nonetheless, the sec- 
ond convention allows us to increase the dimension of 
the system (number of antennas) while maintaining the 
same average received signal to noise ratio (SNR) level 
p ' ct T x , Vi G {1,2}. Moreover, most of the tools from 
random matrix theory used in the asymptotic analysis 
of the achievable data rate of the opportunistic link in 
Sec. IIV-BI require the variance of the entries of channel 
matrices to be normalized by its size. That is the reason 
why the normalized model, i.e., channel transfer matrices 
and power constraints respectively satisfying £[]) and (O, 
was adopted. 

At receiver i G {1,2}, the signal r-j is processed using 
an Ni x Ni matrix Dj to form the N -dimensional 
vector yj = Dji-j. All along this paper, we refer to Dj 
as the post-processing matrix at receiver i. Regarding 
channel knowledge assumptions at the different nodes, 
we assume that the primary terminals (transmitter and 
receiver) have perfect knowledge of the matrix Hn 
while the secondary terminals have perfect knowledge 
of all channel transfer matrices Hy, V(i,j) G {1,2} 2 . 
One might ask whether this setup is highly demanding 
in terms of information assumptions. In fact, there are 



several technical arguments making this setup relatively 
realistic: (a) in some contexts channel reciprocity can be 
exploited to acquire CSI at the transmitters; (b) feedback 
channels are often available in wireless communications 
Ifm . and (c) learning mechanisms lTT2l can be exploited 
to iteratively learn the required CSI. In any case, the 
perfect information assumptions provide us with an 
upper bound on the achievable transmission rate for the 
secondary link. 

III. Interference Alignment Strategy 

In this section, we describe how both links intro- 
duced in Sec. [TT] can simultaneously operate under the 
constraint that no additional interference is generated 
by the opportunistic transmitter on the primary receiver. 
First, we revisit the transmitting scheme implemented by 
the primary system |[T4l . then we present the concept 
of transmit opportunity, and finally we introduce the 
proposed opportunistic IA technique. 

A. Primary Link Performance 

According to our initial assumptions (Sec. J]} the 
primary link must operate at its highest transmission 
rate in the absence of interference. Hence, following the 
results in |[T4l . ifrTI and using our own notation, the 
optimal pre-processing and post-processing schemes for 
the primary link are given by the following theorem. 

Theorem 1: Let Hn = Ue^A/^V^ be a sin- 
gular value decomposition (SVD) of the N\ x M\ 
channel transfer matrix Hn, with ~Uhu an d 
two unitary matrices with dimension N\ x N\ and 
Mi x Mi, respectively, and A# X1 an Ni x Mi matrix 
with main diagonal (A/j^i, . . . , Ajf 11(in i n (jv 1) M; L )) an d 
zeros on its off-diagonal. The primary link achieves 
capacity by choosing Vi = Vir u , Di = U# , 
Pi = diagOi,!, . . . jPi.mJ. where 

Vn€{l,...,M 1 }, Plin =(/3-—?± ) ,(4) 

with, A h h Hii = Ag^Affu 

diag (X h h Hiu1 , . . . , X h h Hii Mi ) and the constant 

(3 (water-level) is set to saturate the power constraint 

(ED- 

Let N = min(iVi, Mi). When implementing its 
capacity-achieving transmission scheme, the primary 
transmitter allocates its transmit power over an equiv- 
alent channel DiHnVi = A^ n which consists of at 
most rank(H^Hn) < N parallel sub-channels with 
non-zero channel gains Xh^Hh ,n> respectively. These 
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non-zero channel gains to which we refer as transmit di- 
mensions, correspond to the non-zero eigenvalues of ma- 
trix H^Hn. The transmit dimension n £ {1, . . . ,M\) 
is said to be used by the primary transmitter if pi jTl > 0. 
Interestingly, (0]) shows that some of the transmit dimen- 
sions can be left unused. Let mi G {1, . . . , Mi} denote 
the number of transmit dimensions used by the primary 
user: 



mi 



Mi 

£ 

n=l 
Mi 

£ 

n=l 



]0,Mipi, max ](Pl,n 



(5) 



,+oo 



(A 



H*H lu n) 



As pi t m ax > 0, the primary link transmits at least over 
dimension n* = are max { A uh h m \ re- 

me{l,...,min(7V 1 ,M 1 )} 1 H ^ Hll ' m > 

gardless of its SNR, and moreover, there exist at most 
A transmit dimensions, thus 



1 < mi < rank(Hf 1 Hn) < A. 



(6) 



In the following section, we show how those unused 
dimensions of the primary system can be seen by the 
secondary system as opportunities to transmit. 

B. Transmit Opportunities 

Once the PA matrix is set up following Th. [T] the 

1/2 1/2 

primary equivalent channel DiHnViP 1 / = A// 11 P 1 / 
is an Ni x M\ diagonal matrix whose main diagonal 
contains mi non-zero entries and A — mi zero entries. 
This equivalent channel transforms the set of mi used 
and Mi — mi unused transmit dimensions into a set of 
mi receive dimensions containing a noisy version of the 
primary signal, and a set of Ai — mi unused receive 
dimensions containing no primary signal. The mi use- 
ful dimensions are called primary reserved dimensions, 
while the remaining Ai — mi dimensions are named 
secondary transmit opportunities (TO). The IA strategy, 
described in Section IIII-CI allows the secondary user to 
exploit these Ai — mi receive dimensions left unused 
by the primary link, while avoiding to interfere with the 
mi receive dimensions used by the primary link. 



Definition 2 (Transmit Opportunities): 



Let 



\ h h # i, . . . \h h flu Mi b e m e eigenvalues of matrix 
H^Hn and (3 be the water-level in (Th. 0). Let mi, 
as defined in (0), be the number of primary reserved 
dimensions. Then the number of transmit opportunities 
S available to the opportunistic terminal is given by 

Mi 

S = N 1 -mi=N 1 



E 1 



,+oo 



"('Vffgtf 11 ,n)- ( 7 ) 



Note that in this definition it is implicitly assumed that 
the number of TOs is constant over a duration equal to 
the channel coherence time. 

Combining © and © yields the bounds on the number 
of transmit opportunities 



Ai - N < S < Ni - 1. 



(8) 



A natural question arises as to whether the number of 
TOs is sufficiently high for the secondary link to achieve 
a significant transmission rate. In order to provide an 
element of response to this question, a method to find 
an approximation of the number of TOs per primary 
transmit antenna, Soo , is proposed in Section IIV-AI In 
any case, as we shall see in the next section, to take 
advantage of the TOs described in this section, a specific 
signal processing scheme is required in the secondary 
link. 



C. Pre-processing Matrix 

In this section, we define the interference alignment 
condition to be met by the secondary transmitter and de- 
termine a pre-processing matrix satisfying this condition. 

Definition 3 (IA condition): Let Hn = 
Uhu Ajfu Vjy be an SVD o/Hn and 

R = a\l Nl +Uf ii H 12 V 2 P 2 VfHf 2 U ifll , (9) 

be the covariance matrix of the co-channel interference 
(CCI) plus noise signal in the primary link. The op- 
portunistic link is said to satisfy the IA condition if its 
opportunistic transmission is such that the primary link 
achieves the transmission rate of the equivalent single- 
user system, which translates mathematically as 



log 2 
log 2 



IjVi+^A^PiAf^ 
IjVi + R- -1 A/in Pi A^ 



(10) 



Our objective is first to find a pre-processing matrix V 2 
that satisfies the IA condition and then, to tune the PA 
matrix P 2 and post-processing matrix D 2 in order to 
maximize the transmission rate for the secondary link. 

Lemma I (Pre-processing matrix V 2 J: Let Hn = 
Uhu Aifn VjJ be an ordered SVD o/H u , with \Jhh 
and Vhh> two unitary matrices of size Ax x Ai and 
Mi x Mi, respectively, and Afj lx an N\ x Mi matrix 
with main diagonal (\hh,i, ■ ■ ■ > AHu,min(iVi,Mi)) and 
zeros on its off-diagonal, such that A^ 1 ^ A^ 2 ^ 
... ^ X 2 tj ■ / Ar ,, r v Let also the Ai x Mo matrix 

^ /in, min(JVi, MiJ 1 z 
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H = H12 have a block structure, 



H 



mi 

Ni — mi 



M 2 
* > 



Hi 

H 2 



(11) 



The IA condition (Def. \3} is satisfied independently of 
the PA matrix P2, when the pre-processing matrix V2 
satisfies the condition: 



HiVo 







miXL 2 ) 



(12) 



where L2 is the dimension of the null space of matrix 
Hi. 

Proof: See Appendix lAl ■ 
Another solution to the IA condition was given in 
ifTTl , namely V2 = H^U^Pi for a given diago- 
nal matrix P\ = diag (pi,i, • • • ,Pi,Afi)> with pi, n = 



(3 I , where (3 is the water-level of the 



primary system (Th. \T} and n G {1, . . . , Mi}. However, 
such a solution is more restrictive than (fT2l since it 
requires Hi 2 to be invertible and does not hold for the 
case when N / Mj, V(i,j) G {1,2} 2 . 
Plugging V2 from (fT2l into © shows that to guarantee 
the IA condition ©, the opportunistic transmitter has to 
avoid interfering with the mi dimensions used by the 
primary transmitter. That is the reason why we refer to 
our technique as OIA: interference from the secondary 
user is made orthogonal to the mi receive dimensions 
used by the primary link. This is achieved by aligning the 
interference from the secondary user with the N\ — mi 
non-used receive dimensions of the primary link. 
From Lemma [T] it appears that the L2 columns of 
matrix V2 have to belong to the null space Ker(Hi) 
of Hi and therefore to the space spanned by the 
dimKer(Hi) = M2 — rank(Hi) last columns of matrix 
Vs , where Hi = U& A a is an SVD of Hnwith 
Ug- and Vg- two unitary matrices of respective sizes 
mi x mi and M2 x M2, and Ag an mi x M2 matrix 
containing the vector (\ 6ll>v . . . , A 5i>inin(mi)M2) ) on its 
main diagonal and zeros on its off-diagonal, such that 
A 2 - > . . . ^ A 2 - i.e., 

Hi ,mm(mi ,M 2 ) 



V 2 G Span (vg nk ^ +1 \ 



(M 2 ) 
Hi 



(13) 



Here, for all i G {1, . . . , M2}, the column vector 
represents the i column of matrix Vg from the left 
to the right. 

In the following, we assume that the L2 columns of the 
matrix V2 form an orthonormal basis of the correspond- 
ing subspace (fT3T ). and thus, V^V2 = Il 2 - Moreover, 



recalling that Hi is of size mi x M2, we would like to 
point out that: 

• When mi < Mi, rank(Hi) < mi and 
dimKer(Hi) > M2 — m\ with equality if and only 
if Hi is full row -rank. This means that there are 
always at least M2 — mi > non-null orthogonal 
vectors in Ker(Hi), and thus, L2 = dimKer(Hi). 
Consequently, V2 can always be chosen to be 
different from the null matrix 0j\/ 2X £ 2 . 
. When, M 2 ^ mi, rank(Hi) < M 2 and 
dimKer(Hi) > 0, with equality if and only if Hi 
is full column-rank. This means that there are non- 
zero vectors in Ker(Hi) if and only if Hi is not 
full column-rank. Consequently, V2 is a non-zero 
matrix if and only if Hi is not full column-rank, 
and again L2 = dimKer(Hi). 
Therefore, the rank of V2 is given by L2 = 
dimKer(Hi) < M2, and it represents the number of 
transmit dimensions on which the secondary transmitter 
can allocate power without affecting the performance of 
the primary user. The following lower bound on L2 holds 

L 2 = dimKer(Hi) = M 2 - rank(Hi) 

> M 2 - min(M 2 , mi) (14) 
= max(0, M2 — mi) 

Note that by processing S2 with V2 the resulting signal 
V2S2 becomes orthogonal to the space spanned by a 
subset of mi rows of the cross-interference channel 
matrix H = H12. This is the main difference 

between the proposed OIA technique and the classical 
zero-forcing beamforming (ZFBF) |fT8l , for which the 
transmit signal must be orthogonal to the whole row 
space of matrix H. In the ZFBF case, the number of 
transmit dimensions, on which the secondary transmitter 
can allocate power without affecting the performance of 
the primary user, is given by L2 ) bf = dimKer(H) = 
M2 — rank(H). Since rank(Hi) < rank(H), we have 
L2,bf < L2. This inequality, along with the observation 
that Ker(H) C Ker(Hi), shows that any opportunity to 
use a secondary transmit dimension provided by ZFBF 
is also provided by OIA, thus OIA outperforms ZFBF. In 
the next section we tackle the problem of optimizing the 
post-processing matrix D2 to maximize the achievable 
transmission rate for the opportunistic transmitter. 

D. Post-processing Matrix 

Once the pre-processing matrix V2 has been adapted 
to perform I A according to (fT3l ). no harmful interfer- 
ence impairs the primary link. However, the secondary 
receiver undergoes the CCI from the primary transmitter. 
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Then, the joint effect of the CCI and noise signals can be 
seen as a colored Gaussian noise with covariance matrix 



Q 



H 2 iV fflI PiVg 1 Hg + a 2 2 l N2 



(15) 



We recall that the opportunistic receiver has full CSI of 
all channel matrices, i.e., Hjj, V(z,j) £ {1,2} 2 . Given 
an input covariance matrix P2, the mutual information 
between the input S2 and the output y 2 = D 2 r2 is 

i? 2 (P 2 ,<r|) = log 2 |l JV2 +D 2 H 22 V 2 P 2 Vf HgDf (D 2 QDf 

^ log 2 |l JV2 +Q-7H 22 V 2 P 2 V 2 H H| f 2 Q-7], (16) 

where equality is achieved by a whitening post- 
processing filter D2 = Q~2 (HI. i.e., the mutual 
information between the transmitted signal S2 and r-2, 
is the same as that between S2 and y 2 = D2r2- Note 
also that expression (fl6l ) is maximized by a zero-mean 
circularly-symmetric complex Gaussian input S2 fffl . 

E. Power Allocation Matrix Optimization 

In this section, we are interested in finding the input 
covariance matrix P2 which maximizes the achievable 
transmission rate for the opportunistic link, R 2 (P 2 ,o 2 ) 
assuming that both matrices V2 and D2 have been set up 
as discussed in Sec. IIII-CI and IIII-DI respectively. More 
specifically, the problem of interest in this section is: 



max log 2 Ijv 2 +Q"3H 22 v 2 P 2 Vf HgQ-s 
P 2 



s.t. 



Trace(V 2 P 2 V 2 H )^p 2iI1 



(17) 



Before solving the optimization problem (OP) in (fTTl ). 
we briefly describe the uniform PA scheme (UPA). 
The UPA policy can be very useful not only to relax 
some information assumptions and decrease computa- 
tional complexity at the transmitter but also because it 
corresponds to the limit of the optimal PA policy in the 
high SNR regime. 

1 ) Uniform Power Allocation: In this case, the oppor- 
tunistic transmitter does not perform any optimization on 
its own transmit power. It rather uniformly spreads its 
total power among the previously identified TOs. Thus, 
the PA matrix P2 is assumed to be of the form 



^2,UPA = llL 2 , 



(18) 



where the constant 7 is chosen to saturate the transmit 
power constraint ©, 



M 2 p 2j] 



M 2 P2,r 



Trace (V 2 Vf) " L 2 



(19) 



2) Optimal Power Allocation: Here, we tackle the 
OP formulated in (fTTl) . For doing so, we assume that 
the columns of matrix V2 are unitary and mutually 

orthogonal. We define the matrix K = Q~2H 2 2V 2 , 
where K is an N 2 x L 2 matrix. Let K = UxA^-Vf 
be an SVD of matrix K, where the matrices JJk and 
V k are unitary matrices with dimensions N 2 x N 2 and 
L 2 x L 2 respectively. The matrix Ak is an N 2 x L 2 matrix 
with at most min(N 2 , L 2 ) non-zero singular values on 
its main diagonal and zeros in its off-diagonal entries. 
The entries in the diagonal of the matrix Ax are denoted 
by Ak,i, ■ ■ ■ , Ax,min(7V 2 ,L 2 )- Finally, the original OP CCZ]> 
can be rewritten as 

argmax log 2 |ljv 2 +AKVgP 2 VKAg| 

Trace(P 2 ) = Trace(vgp 2 v /f ) (20) 

< M 2 p 2 ,max- 



S.t. 



Here, we define the square matrices of dimension L 2 , 

(21) 



P 2 t V f P 2 V K , 



and A K h K = A%A K = diag (X K h K ^, \ k »k,l 2 )- 
Using the new variables P2 and Aj^h K , we can write 
that 



I N2 +A K VgP 2 V K Af 



— \I L2 +A k h k P 2 \ 
L 2 

71=1 



(22) 



where p2, n , with n G {1, . . . , L 2 } are the entries of the 
main diagonal of matrix P 2 . Note that in (1221) equality 
holds if P 2 is a diagonal matrix |[20l . Thus, choosing P2 
to be diagonal maximizes the transmission rate. Hence, 
the OP simplifies to 



max log 2 (l + \K"K,n P2,n) 



P 2 ,1---P 2 ,L 



S.t. 



11 = 1 



(23) 



XI P2,n < M 2 p 2y 



n=l 



The simplified optimization problem d23l has eventually 
a water-filling solution of the form 



Vn €{!,... ,L 2 }, p2,r 



1 



Xkh 



K,n 



, (24) 



where, the water-level (5 2 is determined to saturate the 
power constraints in the optimization problem (l23l . Once 
the matrix P2 (I2TT ) has been obtained using water-filling 
(l24l . we define the optimal PA matrix P2,OPA by 



2,OPA 



diag (p 



2,1 < 



,P2,L 2 ) 



(25) 
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while the left and right hand factors, and V|£, of 
matrix P2 in (|2TT > are included in the pre-processing 
matrix: 

V2,opa = V 2 V K . (26) 

In the next section, we study the achievable transmission 
rates of the opportunistic link. 

IV. Asymptotic Performance of the 
Secondary link 

In this section, the performance of the secondary link 
is analyzed in the regime of large number of antennas, 
which is defined as follows: 

Definition 4 (Regime of Large Numbers of Antennas): 
The regime of large numbers of antennas (RLNA) is 
defined as follows: 

• Vi e {1,2}, Ni -> +00; 

. Vj E {1,2}, Mj -» +00; 



V(i,j) 6 {1,2} 2 , lim 



oqj > is constant. 



M 3 ~>+oc 



Ni 



< +00, and 



A. Asymptotic Number of Transmit Opportunities 

In Sec. Hill two relevant parameters regarding the per- 
formance of the opportunistic system can be identified: 
the number of TOs (S) and the number of transmit 
dimensions to which the secondary user can allocate 
power without affecting the performance of the primary 
user (X2). Indeed, L2 is equivalent to the number of 
idependent symbols the opportunistic system is able to 
simultaneously transmit. In the following, we analyze 
both parameters S and L2 in the RLNA by studying the 
fractions 



L 



2,00 



S 

lim and, 

iVi->+oo M\ 

lim — — . 

Af 2 ->+oo 



(27) 



(28) 



Using (|7]), the fraction can be re-written as follows 

1 



lim — (N ± 

n^+oo Mi 

M 2 ^+oo 



mi, 



1 

an 



mi. 



where, 



mi, 



mi 

hm — — 

N^+oo Mi 
Mi— >+oo 



(29) 



(30) 



As a preliminary step toward determining the expressions 
of and L 2 oo , we first show how to find the asymp- 
totic water-level (3^ in the RLNA, and the expression 



of mi |0O . First, recall from the water-filling solution (O 
and the power constraint (f3]) that 



(31) 



Define the real function q by 
0. 







A 



if A = 0, 
if A > 0, 



which is continuous and bounded on R H 
rewritten as 

Ml 



(32) 
3Tb can be 



1 r°° 

M-Y,*V*&*».») = / f { H f Hl M) dA (33) 

1V±l n=l J -°° 

where fjpijj^ is the probability density function associ- 
ated with the empirical eigenvalue distribution F^ 1 ^^ 
of matrix HfjHn. In the RLNA, the empirical eigen- 
value distribution F^i 1 ^ converges almost surely to the 
deterministic limiting eigenvalue distribution F H n Hll , 
known as the Marcenko-Pastur law |[2ll whose associ- 
ated density is 



f HRHll W = (i-^) s(x) + 



y / (A-a) + (j l -A)+ 



2ttX 



(34) 



where, a =(^i--^=Y and b = (i+-^=Y . Note that the 
Marcenko-Pastur law has a bounded real positive support 
{{0} U [a, b]} and q is continuous and bounded on R+. 
Consequently, in the RLNA, we have the almost sure 
convergence of (l33l) . i.e., 



?(A) 451(A) 



dA 



?(A)/ff gHll (A)<L\ 



Thus, in the RLNA (Def. gj), the water-level fi^ is the 
unique solution |[22l to the equation 



p- 



'-d\- pl 



(35) 



and it does not depend on any specific realization of the 
channel transfer matrix Hn, but only on the maximum 
power pi imax and the receiver noise power a\. 
We can now derive mi qq. From (f5]), we have 



mi, 



1 Mj 
lim — Vl. 

AT^+oo Mi ^ ^,+00 



(AffHffn, 



lim 

JVi — >+oo 
Mi— >+oo 

f6 



1 



(A) C,^) dA 



y/{X-a)(b-X) 



max(a,--J-) 



2ttA 



dA. 



(36) 
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Thus, given the asymptotic number of transmist dimen- 
sions used by the primary link per primary transmit 
antenna mi )0O , we obtain the asymptotic number of 
transmit opportunities per primary transmit antenna 5^ 
by following d27l ), i.e., 



S -J_ 

an 



y/(X-a)(b-\) 



max(a,^i-) 



2ttA 



dA. 



(37) 



From <(8]), the following bounds on Soo hold in the 
RLNA: 



1 



1 < ^oo < 



(38) 



1 

a n j ai2 
Finally, we give the expression of L 2i0 o- Recall that 
L2 = dimKer(Hi) = M2 — rank(Hi). The rank of 
Hi is given by its number of non-zero singular values, 
or equivalently by the number of non-zero eigenvalues of 
matrix Hf Hi. Let Xjjn^ v . . . , ^jj H Hx M 2 denote the 
eigenvalues of matrix H^Hi. We have 

rank(Hi) 



L 



2,00 



1 



lim 

iYi,M 2 — >+oo 



M 2 

M 2 



lim 

N u M a -*+oo 



n=l 
00 



+ oo[(A)/^(A)dA, 
(39) 



where /^f 2 ^ (A) is the probability density function asso- 
ciated with the empirical eigenvalue distribution ■ 

Hi is of size mi x Mi, and the ratio ^ converges in 
the RLNA to 

A M 2 Q 12 



M 2 

ol\ = iim 

JVi,Afi,Af 3 ->oo mi 



< 00. 



anmi oo 



(40) 



Thus, in the RLNA, the empirical eigenvalue distribution 
Fifl 2 } T converges almost surely to the Marcenko-Pastur 
law Ell Fjj H f Ii with associated density 



where c 



1 - — 



'Oil 



5(X) + 
2 



y/(X-c) + (d-\y 



and d 



1 + 



2vrA 
1 



'OLlJ 



(41) 



Using (HJ) in d39]) yields 

"+00 



L 



2,00 



1 



l]o,+oo[(A)/^H^ i (A)dA 



1 {]-cx),o]}(A)/^H^ i (A)dA 
1 

a 1 



(42) 



Thus, given the asymptotic water-level /3oo for the pri- 
mary link, the asymptotic number of TOs per transmit 
antenna is given by the following expression 



L 



2,00 



1 



on 

«12 



-mi. 



(43) 



1 — "li 

«12 



V(A-a)(h-A) 



max(a,--J-) 



2ttA 



dA 



Note that the number (5) of TOs as well as the number 
(L2) of independent symbols that the secondary link 
can simultaneously transmit are basically determined by 
the number of antennas and the SNR of the primary 
system. From (|27T ). it becomes clear that the higher the 
SNR of the primary link, the lower the number of TOs. 
Nonetheless, as we shall see in the numerical examples 
in Sec. |Vj for practical values of SNR there exist a non- 
zero number of TOs the secondary can always exploit. 

B. Asymptotic Transmission Rate of the Opportunistic 
Link 

In this section, we analyze the behavior of the oppor- 
tunistic rate per antenna 



iMP2,<ri)=^log 2 



Iwo+Q^H^VaPaV? H| 



(44) 



in the RLNA. Interestingly, this quantity can be shown 
to converge to a limit, the latter being independent of the 
realization of H22- In the present work, we essentially 
use this limit to conduct a performance analysis of the 
system under investigation but it is important to know 
that it can be further exploited, for instance, to prove 
some properties, or simplify optimization problems |23l . 
A key transform for analyzing quantities associated with 
large systems is the Stieltjes transform, which we define 
in App. |B] By exploiting the Stieltjes transform and 
results from random matrix theory for large systems (See 
App. |B]), it is possible to find the limit of (l44l in the 
RLNA. The corresponding result is as follows. 

Proposition 5 (Asymptotic Transmission Rate): 
Define the matrices 

A 



Mi 

M 2 



M 



A 



A 



H2iV jf / 11 PiVg ii H 2 ^ 

H22V2P2V 2 ¥ H 2 1> 

Mi +M2, 



(45) 
(46) 
(47) 



and consider the system model described in Sec. |77| with 
a primary link using the configuration (Vi, Di, Pi) 
described in Sec. I///-AI and a secondary link with the 
configuration (V 2 , D2, P2) described in Sec. \IITC\ 
I///-DI with P2 any PA matrix independent from the noise 
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level a\. Then, in the RLNA (Def. [?]), under the assump- 
tion that Pi and V2P2V^ have limiting eigenvalue 
distributions Fp 1 and Fy 2 p 2 yH with compact support, 
the transmission rate per antenna of the opportunistic 
link (TX2-RX2) converges almost surely to 

1 r-+oo 

i?2,oc = / G Ml (-z) - G M {-z) dz, (48) 

where, Gm(z) an d Gm 1 (z) are the Stieltjes transforms 
of the limiting eigenvalue distribution of matrices M and 
Mi, respectively. Gm(z) and Gm ± (z) are obtained by 
solving the fixed point equations ( with unique solution 
when z G Br fi2M): 

G M 1 (z) = - ~r, — (49) 



z - g(G Ml (z)) 



and 



G M (z) 



-1 



(50) 



z-g{G M {z))-h{G M {z)Y 

respectively, where the functions g(u) and h(u) are 
defined as follows 



9(u) 



h(u) 




(51) 



(52) 



with the expectations in ri57D and ri52D taken on the 
random variables p\ and P2 with distribution Fp 1 and 
F v 2 P 2 v 2 H > respectively. 

Proof: For the proof, see Appendix O ■ 
The (non-trivial) result in Prop. [5] holds for any power 
allocation matrix P2 independent of a\. In particular, 
the case of the uniform power allocation policy perfectly 
meets this assumption. This also means that it holds for 
the optimum PA policy in the high SNR regime. For 
low and medium SNRs, the authors have noticed that the 
matrix P2,opa is in general not independent of a\. This 
is because P2 is obtained from a water-filling procedure. 
The corresponding technical problem is not trivial and 
is therefore left as an extension of the present work. 

V. Numerical Results 

A. The Number S of Transmit Opportunities 

As shown in (|2"7T ). the number of TOs is a function of 
the number of antennas and the SNR of the primary link. 
In Fig. |2j we plot the number of TOs per transmit antenna 
^oo as a function of the SNR for different number of 
antennas in the receiver and transmitter of the primary 
link. Interestingly, even though the number of TOs is a 
non-increasing function of the SNR, Fig. [2] shows that 



for practical values of the SNR (10 - 20 dBs.) there 
exists a non-zero number of TOs. Note also that the 
number of TOs is an increasing function of the ratio 



(an 



Mi 



). For instance, in the case N\ > M\, i.e., 



an > 1 the secondary transmitters always sees a non- 
zero number of TOs independently of the SNR of the 
primary link, and thus, opportunistic communications are 
always feasible. On the contrary, when an ^ 1, the 
feasibility of opportunistic communications depends on 
the SNR of the primary link. 

Finally, it is important to remark that even though, the 
analysis of the number of TOs has been done in the 
RLNA (Def. H}, the model is also valid for finite number 
of antennas. In Fig. |2l we have also ploted the number 
of TOs observed for a given realization of the channel 
transfer matrix Hn when N\ = 10 and an € {A, 1, 2}. 
Therein, it can be seen how the theretical result from 
(l2~7T ) matches the simulation results. 

B. Comparison between OIA and ZFBF 

We compare our OIA scheme with the zero-forcing 
beamforming (ZFBF) scheme |[T8l . Within this scheme, 
the pre-processing matrix V2, denoted by ^2,zfbf, 
satisfies the condition 



H 12 V 



12 V 2. ZFBF 







N r ,L 2 , 



(53) 



which implies that ZFBF is feasible only in some partic- 
ular cases regarding the rank of matrix H12. For instance, 
when A/2 ^ Ni and H12 is full column rank, the pre- 
processing matrix is the null matrix, i.e., ~Wi,zfbf = 
0m 2 ,l 2 an d thus, no transmission takes place. On the 
contrary, in the case of OIA when M2 ^ N\, it is 
still possible to opportunistically transmit with a non- 
null matrix V2 in two cases as shown in Sec. IIII-CI 

• if mi < M2, 

• or if mi > M2 and Hi is not full column rank. 
Another remark is that when using ZFBF and both 
primary and secondary receivers come close, the oppor- 
tunistic link will observe a significant power reduction 
since both the targeted and nulling directions become 
difficult to distinguish. This power reduction will be less 
significant in the case of OIA since it always holds that 
rank(V2) ^ rank(V2,zFBF) thanks to the existence 
of the additional TOs. Strict equality holds only when 

S = — .As discussed in Sec. IIII-Bl the number 
of TOs (S) is independent of the position of one receiver 
with respect to the other. It rather depends on the channel 
realization Hn and the SNR of the primary link. 
In the following, for the ease of presentation, we consider 
that both primary and secondary devices are equipped 



with the same number of antennas N r 



No 



10 



and Nt = M\ = M2, respectively. In this scenario, we 
consider the cases where N t > N r and N t ^ N r . 

1 ) Case N t > N r : In Fig. we consider the case 
where a ~ §, with N r £ {3,9}. In this case, we 
observe that even for a small number of antennas, the 
OIA technique is superior to the classical ZFBF. More- 
over, the higher the number of antennas, the higher the 
difference between the performance of both techniques. 
An important remark here is that, at high SNR, the 
performance of ZFBF and OIA is almost identical. This 
is basically because at high SNR, the number of TOs 
tends to its lower bound Nt — N r (from ([8])), which 
coincides with the number of spatial directions to which 
ZFBF can avoid intefering. Another remark is that both 
UPA and OPA schemes perform identically at high SNR. 

2) Case N t ^ N r : In this case, the ZFBF solution is 
not feasible and thus, we focus only on the OIA solution. 
In Fig. |4l we plot the transmission rate for the case where 
N r = N t G {3,6,9}. We observe that at high SNR 
for the primary link and small number of antennas, the 
uniform PA performs similarly as the optimal PA. For a 
higher number of antennas and low SNR in the primary 
link, the difference between the uniform and optimal 
PA is significant. To show the impact of the SINR of 
both primary and secondary links on the opportunistic 
transmission rate, we present FigJ5] Therein, it can be 
seen clearly that the transmission rate in the opportunistic 
link is inversely proportional to the SNR level at the 
primary link. This is due to the lack of TOs as stated 
in Sec. IIII-BI For the case when N r < N t with strict 
inequality, an opportunistic transmission takes place only 
if N r — Nt ^ S and Hn is not full column rank. 
Here, the behaviour of the opportunistic transmission rate 
is similar to the case N r = N t with the particularity 
that the opportunistic transmission rate reaches zero at a 
lower SNR level. As in the previous case, this is also a 
consequence of the number of available TOs. 

C. Asymptotic Transmission Rate 

In Fig. [6l we plot both primary and secondary trans- 
mission rates for a given realization of matrices Hjj 
V(z,j) € {1,2} 2 . We also plot the asymptotes obtained 
from Prop. [5] considering UPA in the secondary link and 
the optimal PA of the primary link dU). We observe that 
in both cases the transmission rate converges rapidly to 
the asymptotes even for a small number of antennas. 
This shows that Prop. [5] constitutes a good estimation 
of the achievable transmission rate for the secondary 
link even for finite number of antennas. We use Prop. 
[5] to compare the asymptotic transmission rate of the 
secondary and primary link. The asymptotic transmission 



rate of the primary receiver corresponds to the capacity 
of a single user Nt x N r MIMO link whose asymptotes 
are provided in |[25l . From Fig. [6j it becomes evident 
how the secondary link is able to achieve transmission 
rates of the same order as the primary link depending 
on both its own SNR and that of the primary link. 

VI. Conclusions 

In this paper, we proposed a technique to recycle 
spatial directions left unused by a primary MIMO link, 
so that they can be re-used by secondary links. Interest- 
ingly, the number of spatial directions can be evaluated 
analytically and shown to be sufficiently high to allow a 
secondary system to achieve a significant transmission 
rate. We provided a signal construction technique to 
exploit those spatial resources and a power allocation 
policy which maximizes the opportunistic transmission 
rate. Based on our asymptotical analysis, we show 
that this technique allows a secondary link to achieve 
transmission rates of the same order as those of the 
primary link, depending on their respective SNRs. To 
mention few interesting extensions of this work, we 
recall that our solution concerns only two MIMO links. 
The case where there exists several opportunistic devices 
and/or several primary devices remains to be studied in 
details. More importantly, some information assumptions 
could be relaxed to make the proposed approach more 
practical. This remark concerns CSI assumptions but 
also behavioral assumptions. Indeed, it was assumed 
that the precoding scheme used by the primary trans- 
mitter is capacity-achieving, which allows the secondary 
transmitter to predict how the secondary transmitter is 
going to exploit its spatial resources. This behavioral 
assumption could be relaxed but some spatial sensing 
mechanisms should be designed to know which spatial 
modes are effectively used by the secondary transmitter, 
which could be an interesting extension of the proposed 
scheme. 

Appendix A 
Proof of Lemma Q] 

Here, we prove Lemma [7] which states that: if a matrix 
V2 satisfies the condition H1V2 = ^{n 1 -s)xL 2 men ^ 
meets the I A condition ©. 

Proof: Let Hn = Uh^Ajj^VJj be a sorted 
SVD of matrix Hn, with TJhu an d Vjy 1I5 two uni- 
tary matrices of sizes Nt x N\ and M\ x M\, re- 
spectively, and A.Hn an -^1 x Mi matrix with main 
diagonal (\ Hll ,i, • • • , A frilimilI ( J vk I JMi)) and zeros on its 
off-diagonal, such that 1 ^ A^ 2 ^ • • • ^ 
^Hn min(JVi Mi) - Given that the singular values of the 



1 1 



matrix Hu are sorted, we can write matrix Afl^PiA^ 
as a block matrix, 



LXtJVi-mi) 



0, 



'(ifl-mi)Xmi 0(N 1 -,T 11 )x(Ni-m 1 ) 

where the diagonal matrix \I/ of size mi x mi is * 



(54) 



Now let us split the interference-plus- 



split 
matrix © as: 



plus-noise covanance 



N 1 -m 1 



R= 



'"l 



Ri+(T 2 I mi 



R 2 

R; 



-cr 2 I Nl - 



(55) 

where (Ri + crfl mi ) and (R 3 + afl Nl _ mi ) are invert- 
ible Hermitian matrices, and matrices Ri, R2 and R3 
are defined from © and (TTTb as 



Ri 4 H!V 2 P 2 VfHf, 
R2 — HiV2P2V 2 ^H 2 ¥ , 
R3 = HoV2P2V?H^. 



(56) 
(57) 

i 2 v 2 ±-2 v 2 n 2 • ( 58 ) 

Now, by plugging expressions (l54l and (I55T ) in (TTOb . the 
IA condition can be rewritten as follows: 



log 2 |<r 2 I mi +*| 

— 10g 2 |Rl+(TfI; 

l°g 2 ' 



log 2 |o- 2 I Nl |=log 2 |Ri+o-?I mi +*| 
n - R ¥ (R.i+^?i mi ) _1 R. 2 | 



Rs + 'l^-m! -R| f (Rl+ CT flm 1 +*) R-2 



(59) 



Note that there exists several choices for the submatrices 
Ri, R2, and R3 allowing the equality in d59l to be met. 
We see that a possible choice in order to meet the IA 
condition is Ri = 0, R2 = 0, independently of the 
matrix R3. Thus, from d56l ) and d57l ) we have Ri = 
and R2 = by imposing the condition H1V2 = 
miX L 2 , for any given PA matrix P2, which concludes 
the proof. ■ 

Appendix B 
Definitions 

In this appendix, we present useful definitions and 
previous results used in the proofs of Appendix O 

Definition 6: Let X be an n x n random matrix 

(ri) 

with empirical eigenvalue distribution function F x . 
We define the following transforms associated with the 
distribution for z G C+ = {z G C : Im(z) > 0}: 

A f°° 

Stieltjes transform: G x (z) = I j^dF^\t), (60) 



r x (z) 



A 



A 



/oo 
j^dF^ W ,(61) 
-00 



S-transform: S x (z) = ■^T x 1 ( 



(62) 



where the function T x (z) is the reciprocal function of 
T x (z), i.e., 

T x 1 (T x (z)) = T x (T x 1 (z)) = z. (63) 

From (l60l and d6TT >. we obtain the following relationship 
between the function T x {z) (named T-transform in 
||26l ) and the Stieltjes transform G x (z), 



1 



T x {z) = -l--G x 



(64) 



Appendix C 
Proof of Proposition [5] 

In this appendix, we provide a proof of Prop. [5] on the 
asymptotic expression of the opportunistic transmission 
rate per antenna, defined by 



£2,00 (P 2 ,0 = lim 

V(i,i)e{l,2} 2 , N^M^oo 
V(ij)e{l,2} 2 , ^->a„<oc 



J R 2 (P 2 ,CT 2 



First, we list the steps of the proof and then we present 
a detailed development for each of them: 

1) Step 1: Express dfl 2 .oc(P 2 ,o- 2 ) as function of the 
Stieltjes transforms Gmi(z) and Gm{z), 

2) Step 2: Obtain G Ml (z), 

3) Step 3: Obtain G M (z), 

4) Step 4: Integrate 8 ^ 2 - ^ 2 ' g ' } to obtain 

-R2,oo(P2, _ 

Step 1: Express gga^gnd) as a function of the 
Stieltjes transforms Gm x {z) and Gm{z)- 
Using (fT6l ) and ([TBI ), the opportunistic rate per receive 
antenna R? can be re-written as follows 



J f?2(P 2 , ' 2 2 ) 



^log 2 Ijv 2 +Q-JH 22 V 2 P 2 VfH|Q- 



N 2 

A 



1JV 2 

log 2 [ <j\ I N2 +M 1 +M 2 1 - ^ log 2 



with Mi H 2 iV fl - 1I PiV^ ii H|i, M 2 

H22V 2 P 2 Vf ng, and M = m"+ M 2 . Matrices 
M and Mi are Hermitian Gramian matrices with 



in 



(65) 

-Mi|, 

A 



U M A M U 



H 



eigenvalue decomposition M = 
and Mi = U^A^U]^ , respectively. 
Matrix Um and Vmi are A2 x N2 unitary 
matrices, and A M = diag(AAf,i, . . . , Xm,N 2 ) and 
A Ml = diag(A A // li i,... ,Aa/ 1iA t 2 ) are square diagonal 
matrices containing the eigenvalues of the matrices M 
and Mi in decreasing order. Expression (l65l) can be 
written as 

N 2 

R 2 (?2,vl) = ^J^loga^+AM.O-kfeCoS+AMi.O (66) 



= 1 



(A+<r 2 )dF<f 2 » (A)— log 2 {X+al )dFl" 2 > ( A) 



(N 2 ) , 



log 2 (A+ CT2 2 )dF M (A)- / log 2 (A+a 2 )dF Ml (A) 
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where F^ 2 ' and F^ 2 ^ are respectively the empirical 
eigenvalue distributions of matrices M and Mi of size 
N2, that converge almost surely to the asymptotic eigen- 
value distributions Fm and Fm ± , respectively. Fm and 
Fm 1 have a compact support. Indeed the empirical eigen- 
value distribution of Wishart matrices H^-H^ converges 
almost surely to the compactly supported Marcenko- 
Pastur law, and by assumption, matrices VjPjV,^, i G 
{1,2} have a limit eigenvalue distribution with a com- 
pact support. Then by Lemma 5 in E71 , the asymptotic 
eigenvalue distribution of Mi and M2 have a compact 
support. The logarithm function being continuous, it is 
bounded on the compact supports of the asymptotic 
eigenvalue distributions of Mi and M, therefore, the 
almost sure convergence in d66l ) could be obtained by 
using the bounded convergence theorem ll28l . 
From d66l ), the derivative of the asymptotic rate 
^2,00 (P2 5 f 2 ) with respect to the noise power a\ can 
be written as 



r ii2,oo(P2,0i) 



■dF Ml (A) 



(67) 



where Gm {z) and Gm x {z) are the Stieltjes transforms 
of the asymptotic eigenvalue distributions Fm and Fm 1 , 
respectively. 

Step 2: Obtain G Ml {z) 

Matrix Mi can be written as 

Mi = v^EbiV^^Vf^Hliv^i". (68) 

"21 

The entries of the JV 2 x M\ matrix y / a 2 7H 2 i are zero- 
mean i.i.d. complex Gaussian with variance = 
thus y / a2iH2i is bi-unitarily invariant. Matrix ~Vhu is 
unitary, consequently ^/a^iH^iVfl-^ has the same distri- 
bution as y / aiiH2i, in particular its entries are i.i.d. with 
mean zero and variance From (01), is diagonal, 
and by assumption it has a limit eigenvalue distribution 
Fpj . Thus we can apply Theorem 1.1 in 11241 to Mi, in 

the particular case where A = 0tv 2 to obtain the Stieltjes 
transform of the asymptotic eigenvalue distribution of 
matrix Mi 



G Ml (z) 



Go N \z-a 21 ___ Q , 21 / fJi ( Q , 21 A)dA 



G (z-g{G M ^{z))), 



(69) 



where the function g(u) is defined by 

/oo 
TT ^ r / Pl (t)dt = E 
-OO "21 



H L —tu 



where the random variable t follows the c.d.f. Fp 1 . 
The square null matrix has an asymptotic eigenvalue 
distribution Fq{\) = //(A). Thus, its Stieltjes transform 
is 

Go(z) = I" — *— 5(A)dA = --. (70) 

./ — no A — Z Z 



Then, using expressions 



and (1701) . we obtain 
-1 



(71) 



z - g(G Ml (z))' 

Expression (1711 ) is a fixed-point equation with unique 

solution when z G IRT IT241 . 

Step 3: Obtain G M (z) Recall that 

M 4 H 22 V 2 P 2 Vf Hg + H2iV Hll PiVg ii H^ 1 (72) 

To obtain the Stieltjes transform Gm, we apply Theo- 
rem 1.1 in 11241 as in Step 2: 



G M (z) = GmAz-9(Gm(z))). 



(73) 



To obtain the Stieltjes transform Gm 2 °f me asymptotic 
eigenvalue distribution function of the matrix M2 = 
H22V2P2V^H2^, we first express its ^-transform as 

Sm 2 {z) = S H22V2 p 2V h h h(z) 

= S/^H 22 V2^-V 2 H H 2 H 2 ^/a 2 -2~( Z )> 



and by Lemma 1 in l[27\l : 



ij\/ 2 



. Z + «22 J S V^2 2 » 22 H 22 V^ 22 V 2 -Si- V,# ( Q2 2 ^' 



and by Theorem 1 in l\29\ 



Sm 2 (z) 



. z + » 22 J S ^22hH H22 v^22 (, c 22 J S v 2 V. 1 



8+1 



i + Q 22 / \^l+0 22 ^^ 



(74) 



The S-transforms Sm 2 (z) and 5y 2 p 2 yH (^) in expres- 
sion d74b can be written as functions of their T- 
transforms: 



(75) 



s M2 {z) = ±±*T^(>),from ()62 



'm 2 v,i — - ~ j -m 2 



S..Avh(±) = from (1621) 



2+£ T -l 

v 9 -S-v» v Q 2 2 y 

=22 



Then, plugging (1751 ) and d76l ) into (1741 yields 



l + JV^ ( a22 ) 



(76) 



(77) 
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Now, using the relation ( f64T > between the T-transform 
and the Stieltjes transform, we write 

GmM = (^pj (?M 2 Q) + l) , (78) 



and from (|73T ). we obtain 



G M (z) 



z-g(G M (z))) \ T M 2 [ z - g (G M (z))) + 1 ) -( 79 ) 

We handle ( p79T > to obtain Gm(z) as a function of 

^V 2 P2Vf( z ) : 



! U-!(G„(«)) 



= -1-(2- 9 (Gm(z)))Gm(2) (80) 

^(Tyy = T^ 2 {-l-(z-g(G M (z)))G M {z)) 



- 1 



-s(G M (z)) ('-j(U|kW))GikW 

T -l ( l + (z~g(G M (z)))G M (z) ^ 

l + (z-g(G M (z)))G M (z) \ 



■Gjif(z) 



x 



2_ V H 



V 2 — <i- V, 
"22 ' 

l + ( i-jIGmMI ) G M (z) 



'22 



-s(G M (z)), 



G«{z) 



From the definition of the T-transform (1611 ), it follows 
that 



l+a 2 2T v2 ^ 2 _ vH (-G M (z)) ). 



Q2 2 T^ p 2 vH (-Gm(z)) 



a-22 



It 



G M (z)\ 



+ G M (z)X v 2 

= / 7X")^ a22 ^ 2f - 2 v^(^A)dA 

= r - g mm\ f VPvH (t)dt (81) 

J l + G M (z)^L- J V 2 P 2 V» W V / 

Using (HB in ([80]>, we have 

G "M=( ^-»i%w) ) (i-GmW MGmW) ) (82) 
with the function /i(u) defined as follows 

P2 



h(u) 



I Y-^- t dF v 2 p 2 v-(t)=^ 



1 + 77~P2 U 



where the random variable p2 follows the distribution 

f v 2 p 2 v 2 h - 

Factorizing Gm{z) in d82l finally yields 

G M {z) 



z-g(G M (z))-h(G M (z)) (83) 

Expression (l83l) is a fixed point equation with unique 
solution when z € R- ll24l . 

Step 4: Integrate dR2 f^ to obtain i? 2 (P 2 ,<7i) in 
the RLNA. 

From d67l ), we have that 



r J R2,o (P2,<T 2 2 )= I ^(GM(-<Tl)-GM 1 (- CT l)). 



(84) 



Moreover, it is know that if a\ — > od no reliable 
communication is possible and thus, -fi^.oo = 0. Hence, 



the asymptotic rate of the opportunistic link can be 
obtained by integrating expression (l84l 



-1 f°° 

R2 '°° = ir72 / 2 (yGhI (~ Z ^ ~ ° Ml (~ Z ^ dz ' (85) 
which ends the proof. 



References 

[1] S. M. Perlaza, M. Debbah, S. Lasaulce, and J.-M. Chaufray, 
"Opportunistic interference alignment in MIMO interference 
channels," in Proc. IEEE 19th Intl. Symp. on Personal, Indoor 
and Mobile Radio Communications (P1MRC), Cannes, France. 
Sept. 2008. 

[2] S. Haykin, "Cognitive radio: brain-empowered wireless com- 
munications," IEEE J. Select. Areas Commun., vol. 23, no. 2, 
pp. 201-220, Feb. 2005. 

[3] M. Haddad, A. Hayar, and M. Debbah, "Spectral efficiency of 
spectrum-pooling systems," IET Commun., vol. 2, no. 6, pp. 
733-741, July 2008. 

[4] V. Cadambe and S. Jafar, "Interference alignment and degrees 
of freedom of the A'-user interference channel," IEEE Trans. 
Inform. Theory, vol. 54, no. 8, pp. 3425-3441, Aug. 2008. 

[5] S. Jafar and M. Fakhereddin, "Degrees of freedom for the 
MIMO interference channel," IEEE Trans. Inform. Theory, 
vol. 53, no. 7, pp. 2637-2642, July 2007. 

[6] M. Maddah-Ali, A. Motahari, and A. Khandani, "UW-ECE- 
2006-12 - communication over X channel: Signalling and 
multiplexing gain," University of Waterloo, Tech. Rep., 2006. 

[7] , "Communication over MIMO X channels: Interference 

alignment, decomposition, and performance analysis," IEEE 
Trans. Inform. Theory, vol. 54, no. 8, pp. 3457-3470, Aug. 
2008. 

[8] H. Weingarten, S. Shamai, and G. Kramer, "On the compound 
MIMO broadcast channel," in Proc. Annual Information Theory 
and Applications Workshop, San Diego, USA. Feb. 2007. 

[9] V. Cadambe, S. Jafar, and S. Shamai, "Interference alignment 
on the deterministic channel and application to fully connected 
AWGN interference networks," in Proc. IEEE Information 
Theoty Workshop (ITW), Porto, Portugal. May. 2008. 
[10] S. W. Peters and R. W. Heath, "Interference alignment via al- 
ternating minimization," in Proc. IEEE Intl. Conf. on Acoustics, 
Speech and Signal Processing (ICASSP), Taipei, Taiwan. Apr. 
2009. 

[11] J. Thukral and H. Bolcskei, "Interference alignment with limited 
feedback," in Proc. IEEE Intl. Symp. on Information Theory 
(ISIT), Seoul, Korea. Jun. 2009. 

[12] K. Gomadam, V. Cadambe, and S. Jafar, "Approaching the 
capacity of wireless networks through distributed interference 
alignment," in Proc. IEEE Global Telecommunications Confer- 
ence (GLOBECOM), New Orleans, USA. Dec. 2008. 

[13] R. Tresch, M. Guillaud, and E. Riegler, "On the achievability 
of interference alignment in the K-user constant MIMO inter- 
ference channel," in Proc. IEEE Workshop on Statistical Signal 
Processing (SSP), Cardiff, Wales, UK. Sept. 2009. 

[14] E. Telatar, "Capacity of multi-antenna Gaussian channels," 
European Transactions on Telecommunications, vol. 10, no. 6, 
pp. 585-596, 1999. 

[15] S. Vishwanath and S. Jafar, "On the capacity of vector Gaussian 
interference channels," in Proc. IEEE Inform. Theoty Workshop 
(ITW), San Antonio, USA. Oct. 2004. 

[16] X. Shang, B. Chen, and M. Gans, "On the achievable sum rate 
for MIMO interference channels," IEEE Trans. Inform. Theory, 
vol. 52, no. 9, pp. 4313^1320, Sept. 2006. 



14 



[17] 
[18] 

[19] 

[20] 
[21] 

[22] 
[23] 

[24] 

[25] 

[26] 
[27] 

[28] 
[29] 



E. Telatar, "Capacity of multi-antenna Gaussian channels," Bell 
Labs., Tech. Rep., 1995. 

A. Paulraj, R. Nabar, and D. Gore, Introduction to Space-Time 
Wireless Communications. Cambridge, U.K: Cambridge Univ. 
Press, 2003. 

F. Neeser and J. Massey, "Proper complex random processes 
with applications to information theory," IEEE Trans. Inform. 
Theory, vol. 39, no. 4, pp. 1293-1302, Jul 1993. 

A. W. Marshall and I. Olkin, Inequalities: Theory of Majoriza- 

tion and Its Applications. Academic Press, 1979. 

V. A. Marcenko and L. A. Pastur, "Distribution of eigenvalues 

for some sets of random matrices," Mathematics of the USSR- 

Sbornik, vol. 1, no. 4, pp. 457^183, 1967. 

C. E. Chuah, D. N. C. Tse, J. M. Kahn, and R. A. Valenzuela, 

"Capacity scaling in MIMO wireless systems under correlated 

fading," IEEE Trans. Inform. Theory, vol. 48, pp. 637-650, 

2002. 

J. Dumont, W. Hachem, S. Lasaulce, P. Loubaton, and 
J. Najim, "On the capacity achieving covariance matrix of 
Rician MIMO channels: An asymptotic approach," to appear 
in IEEE Trans, on Inform. Theory: Revised on 2009. [Online]. 
Available: http://arxiv.org/abs/0710.4051 

J. W. Silverstein and Z. D. Bai, "On the empirical distribution 
of eigenvalues of a class of large dimensional random matrices," 
J. Multivar. Anal., vol. 54, no. 2, pp. 175-192, 1995. 

G. J. Foschini, "Layered space-time architecture for wireless 
communication in a fading environment when using multi- 
element antennas," Bell Labs, Tech. Rep., 1996. 

Tulino and S. Verdu, Random Matrix Theory and Wireless 

Communications. Now Publishers Inc., 2004. 

N. Fawaz, K. Zarifi, M. Debbah, and D. Gesbert, "Asymptotic 

capacity and optimal precoding strategy of multi-level 

precode & forward in correlated channels," Submitted to 

IEEE Trans. Inform. Theory, 2008. [Online]. Available: 

|http://arxiv.org7 abs/0906.0037l 

R. G. Bartle, The Elements of Integration and Lebesgue Mea- 
sure. Wiley-Interscience, 1995. 

R. Miiller, "On the asymptotic eigenvalue distribution of con- 
catenated vector-valued fading channels," IEEE Trans. Inform. 
Theory, vol. 48, no. 7, pp. 2086-2091, July 2002. 



Primary System 




Rx 2 



H 



H 



22 



21 X 



Secondary System 




Figure 1. Two-user MIMO interference channel. 
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Figure 2. Fraction of transmit opportunities in the RLNA (Def. 
ED, i.e., Soo © as function of the SNR = ^f^- and a u = ^, 
Simulation results are obtained by using one realization of the matrix 
Hn when Ni - 10. 
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Figure 3. Transmission rate of the opportunistic link obtained by 
Monte Carlo simulations as a function of the SNRi = SNR2 when 
IA and ZFBF are implemented. The number of antennas satisfy a = 
^ w |, with Mi = Mi = N t and Ni = N 2 = N r G {3, 9} and 
SNRi = 2*1^5*, for all i £ {1, 2}. 
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Figure 4. Transmission rate of the opportunistic link obtained by 
Monte Carlo simulations as a function of the SNRi = SNR2. The 
number of antennas satisfy Mi = M2 = N t and Ni = N2 = N r , 
with Nt = N r , and N r € {3,6,9} and SNRi = for all 

%€ {1,2}. 




Figure 5. Transmission rate of the opportunistic link obtained by 
Monte Carlo simulations as a function of the SNRi = p '^ ax , with 
i £ {1,2}. The number of antennas satisfy Mi = M2 = Nt and 

iVi = N 2 = N r , with iV r = N t = 4. 
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Figure 6. Asymptotic transmission rates per antenna of the oppor- 
tunistic link as a function of the number of antennas when N r = Nt 
using uniform PA at different SNR levels SNRi = p ^*| ax . Simulation 
results are obtained using one channel realization for matrices 
V(i, j) £ {1,2} 2 and theoretical results using Prop. 



